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Abstract. Developed magnetohydrodynamic turbulence near two dimensions 
d up to three dimensions has been investigated by means of renormalization 
group approach and double expansion regularization. A modification of standard 
minimal subtraction scheme has been used to analyze the stability of the 
Kolmogorov scaling regime which is governed by the renormalization group fixed 
point. The exact analytical expressions have been obtained for the fixed points. 
The continuation of the universal value of the inverse Prandtl number u = 1.562 
determined at d = 2 up to d = 3 restores the value of u = 1.393 which is known 
in the kinetic fixed point from usual e-expansion. The magnetic stable fixed point 
has been calculated and its stability region has been also examined. This point 
losses stability: (1) below critical value of dimension d c = 2.36 (independently 
on the a-parameter of a magnetic forcing) and, (2) below the value of a c = 0.146 
(independently on the dimension). 
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1. Introduction 

The renormalization group (RG) methods have been widely used to the analysis of 
fully developed hydrodynamic (HD) turbulence beginning from pioneering papers [1, 2] 
based on [3, 4]. It gives possibility to reply upon some principal questions, e.g., on 
the fundamental description of the infrared (IR) scale invariance, as well as it is useful 
for calculation of many quantities, e.g., critical dimensions of the fields and their 
gradients, viscosity, etc. (see, e.g. [5, 6, 7, 8]). 

Then many authors begin to use Wilson's scheme or some adequate 
generalized renormalization scheme to study of HD turbulence [9] as well as of 
magnetohydrodynamic (MHD) turbulence [10, 11]. This time Vasiliev's team have 
used functional formulation of the field-theoretic RG [12, 13] to legalize of the 
Kolmogorov scaling regime of HD turbulence [14, 15]. They consider (as used in 
present paper) the functional quantum field RG approach [5] rather then Wilson's RG 
technique [16]. It assigns a field action to the stochastic problem and makes possible 
to use elegant and very well developed RG procedure in quantum field theory to 
investigate infrared asymptotic regimes of a stochastic system. Then this RG method 
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has been applied in MHD turbulence [17, 18]. Note here that this functional RG 
method allow a straightforward extension of the pertubative calculation to an higher 
order loops without a principal difficulty (see [19, 20], for example). 

Considerable effort had been devoted to application of adequate field-theoretical 
methods in the MHD turbulence (see recent review of Verma [21], for example). 
Authors in [10, 11] have used the 'classical' Yakhot-Orszag scheme [9]. In last years 
Verma [22, 23] performed detailed RG calculation of MDH turbulence using McComb's 
alternative field-theoretic RG procedure [7] and they reached notable progress in 
calculation of some renormalized parameters of MHD turbulence. Here we will not 
present full discussion of all methods used in the full developed turbulence theory 
such as calculation of Alfven ratio, magnetic resistivity [23] or a problem of magnetic 
dynamo in helicit MHD [18] because it goes out of the frame of present paper (but 
see some remarks and discussion in section 4 and section 5). 

Present paper deals with an investigation of the existence and range of stability 
of the 'magnetic' scaling regime (i.e. the magnetic fixed point for zero inverse Prandtl 
number, see below) in the non-helical d-dimensional MHD turbulence. The existence 
of two different anomalous scaling regimes in three dimensions, which are known as 
kinetic and magnetic ones, was established in the pioneering papers [10, 17]. These 
two points correspond to two IR stable fixed points of the RG. On the other hand, 
it was also supposed that in two dimensions the magnetic fixed point does not exist 
as a result of nonexistence of the IR stable magnetic fixed point. But the conclusions 
about two dimensional fixed points cannot be consider without doubts in these papers 
due to the problems with renormalization in two dimensions which were not taken into 
account [24] (see also [5]). In [25, 26] two dimensional case was studied too but again 
with shortcomings, therefore their results cannot be considered completely conclusive. 
Within the our field theoretic RG approach the problem is related to the existence of 
additional divergences which arise in two dimensions. 

The first correct treatment of the two dimensional case of the stochastically forced 
MHD equations with the proper account of these additional divergences was done 
in [27]. It was accomplished within a two-parameter expansion (double expansion) 
of scaling exponents and scaling functions [24] where, besides the parameter which 
characterize the deviation of the exponent of the powerlike correlation function of 
random forcing from its critical value, the additional parameter of the deviation of 
the spatial dimension was introduced. The using of this double expansion method has 
allowed them to confirm the basic conclusions of the previous works [10, 17], namely, 
the nonexistence of the magnetic scaling regime near two dimensions. 

The authors of the paper [27] also tried to restore the stability of the magnetic 
fixed point when moving from two dimensions in direction of three dimensions. This 
possibility was achieved by using of the special choice of finite renormalization which 
allowed them to keep track of the effect of the additional divergences near two 
dimensions. Technically, it was done by introducing of another uniform UV cutoff 
in all propagators which does not affect the large-scale properties of the model. This 
setup is similar to that of Polchinski [28]. As a result, the borderline dimension between 
stable and unstable magnetic fixed point was found and it leads to the possibility of 
the uniform description of two and three dimensional cases of stochastic MHD. 

Another possibility how to solve the problem of the additional divergences in 
two dimensions together with the problem of restoration of the stability of the 
corresponding fixed point when going from a two dimensional system to a three 
dimensional one was proposed in [29]. They suggest to apply a modified minimal 
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substraction (MS) scheme in which the d-dependence of the tensor structures of the 
UV divergent parts of the corresponding diagrams are kept. It was successfully used in 
the fully developed Navier-Stokes turbulence with weak uniaxial anisotropy to restore 
the stability of the Kolmogorov scaling regime which is unstable in two dimensions 
and stable in three dimensions. 

In what follows, we shall apply the double expansion method together with 
modified MS scheme introduced in [29] to the stochastic MHD equations. Our aim is to 
investigate if it is possible to describe correctly and uniformly the two dimensional and 
the three dimensional systems and to compare our results with that of [27] where the 
different method was used (see above). Thus, we carry out an analysis of the randomly 
forced MHD equations with the proper account of the additional UV divergences which 
are appeared in d — 2. We apply a modified minimal subtraction scheme based on 
the fact that the tensor structure of counter-terms is left generally d-dependent in the 
calculations of divergent parts of Green's functions. It will be shown that it allows us 
to investigate behavior of the system under continual transition to d = 3 beginning 
from d — 2. We have also confirmed the earlier conclusions made in [10, 17, 22] that 
near two dimensions a scaling regime driven by the velocity fluctuations may exist, 
but no magnetically driven scaling regime can occur. We have also investigated the 
long-range asymptotic behavior of the model in the double expansion framework and 
found, in particular, that in this case thermal fluctuations of the magnetic scaling 
regime may occur and that the value of the borderline dimension is significantly lower 
(d c = 2.36) than in the e expansion [10] (d c — 2.85) and rather lower than in the 
'modified' double expansion introduced in [27] (d c = 2.46) but it is rather higher 
then value (d c = 2.2) calculated in the frame of the McComb's renormalization [21]. 
The discrepancy between the value of inverse Prandtl number u which corresponds 
to nontrivial stable fixed point of the RG in the three dimensions, which has been 
obtained in the double expansion scheme in earlier paper [30] and that obtained by 
the usual e-expansion scheme [10, 17] and also that obtained by Verma [22, 23] by 
McComb's procedure, was one more reason of the present analysis. Here we show 
that the continuous transition from the universal value of the inverse Prandtl number 
u = 1.562 determined at d = 2 restores the value of u = 1.393 at d = 3 which is known 
from usual e-expansion. 

The paper is organized as follows: In section 2 the functional field theoretic 
formulation of the model is present in detail. In section 3 the renormalization of the 
model is discussed. In section 4 detailed analysis of the possible scaling regimes is 
done. In section 5 conclusions and discussion of the results are given. 

2. Functional formulation of double expansion model 

In the present paper we study the universal statistical features of the model of 
stochastic MHD which is described by the system of equations for the fluctuating 
velocity field of an incompressible conducting fluid v(x), x = (x,t), V-v = and 
the magnetic induction B = (pfi) 1 ^ 2 b(x) (where p is density of the fluid and \i is its 
permeability) [10, 17, 31]: 



with the incompressibility conditions V /" = and V-/ 6 = and the field b 
is suppose to be solenoidal too, V-fo = 0. The statistics of v, b is completely 



d t v + {v-V)v - (6-V)6 - vqS7 2 v = f v , 
d t b + (v-W)b - (b-V)v - v u V 2 b = f , 



(1) 
(2) 
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determined by both the non-linear equations (1,2) and the statistics of the external 
inter-correlated large-scale random forces f b . The dissipation is controlled by the 
parameter of kinematic viscosity v Q , and u denotes inverse Prandtl number (hereafter 
all parameters with a subscript denote bare parameters of unrenormalized theory; 
see below). Note here that the term (b-V)b expresses the transverse part of Lorentzian 
force and it can be omitted in the case of magnetic field treated as a passive admixture. 

As usual [10, 17], statistical properties of the Gaussian forcing with zero mean 
values ((/")= 0, ( f b ) = 0) are determined by relations: 

(//(1)/J(2)) = 5{r)u vl J -0- d V js (k)j k -* [g vW k 2 - 2S - 2 * + g v20 k 2 ] , (3) 

(/j-(l)/ s 6 (2)) =5{r)ulvl J ^P js (fe)e ife - [.g fcl0 fc 2 "" + .g fc20 fc 2 ] , (4) 

where the argument 1 = x\, t = t\ — <2, x = x x — x 2 , Pj S (k) = Sj S — kjk s /k 2 , the 
parameter e determines the powerlike falloff of the long-range forcing correlations, 
and the parameter 5 = (d — 2)/2 describe the deviation from spatial dimension 
d = 2. The free parameter a controls the power form of magnetic forcing. Note 
that parameters e = 2, a = 1 are the natural "physical" values in our "massless" 
power-law energy injection. The introduction of the local correlations (proportional 
to the new couplings g V 2o, and g&2o) which are described by the analytic in k 2 terms 
in the correlation functions (3), and (4) is related to the existence of additional 
divergences of this structure (see below in the text) in the two dimensional model 
which cannot be removed by corresponding nonlocal terms [24, 32, 33]. At the same 
time, the localness of the counterterms is the fundamental feature of a model to be 
multiplicatively renormalizable [34, 13]. For example, it was not taken into account 
in the analysis of the model in [10, 17]. 

Using the well-known Martin-Siggia-Rose formalism [3, 4], one can transform the 
stochastic problem (l)-(2) with correlators (3), and (4) into the field theoretic model 
of the doubled set of fields $ = {v, b, v', b'} with the following action functional 

S=\j dx,dx 2 {^(1) <#(l)/ a "(2)>o«i(2) + ^.(l)(/j , (l)/ s 6 (2)) 6U2)} + 
+ J dx v'- (-d t v + vq V 2 v - (v-V)v - (b-V)b) 

+ J dx b'- (-d t b + u ^ V 2 6 + (6-V)« - (v-V)b) , (5) 

where v', and b' are independent of v, and b auxiliary incompressible fields, which we 
have to introduce when transforming the stochastic problem into a functional form. 

The dimensional constants g v io, g&io, g V 20, and gb20, which control the amount of 
randomly injected energy given by (3), (4), play the role of the coupling constants 
in the perturbative expansion. For the convenience of further calculations the factors 
VqUq and v^Uq including the "bare" (molecular) viscosity vq and the "bare" (molecular 
or microscopic) magnetic inverse Prandtl number u have been extracted. As was 
mentioned already the bare (non-renormalized) quantities are denoted by subscript 
"0". 

The most important measurable quantities in the study of a fully developed 
turbulence and related problems are considered to be the statistical objects represented 
by correlation and response functions (Green functions) of the fields. Standardly, 
the formulation through the action functional (5) replaces the statistical averages of 
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random quantities in the stochastic problem (l)-(4) with equivalent functional averages 
with weight expS(<I>). Generating functionals of total Green functions G(A) and 
connected Green functions W(A) are then defined by the functional integral 

G(A) = e w ^ = Jv<S> e *(*)+^, (6) 

where A(x) = {A v , A b , A v , A b } represents a set of arbitrary sources for the set of 
fields <&, £>$ = VOVO'VvVv' denotes the measure of functional integration, and linear 
form A<S> is defined as 

A<$>= ( dx[A v (x) ■ v(x) + A h (x) ■ b(x) + A v '(x) ■ v'(x) + A b '(x) • b'{x)]. (7) 



The functional formulation gives the possibility of using the field theoretic 
methods, including the RG technique to solve the problem. By means of the RG 
approach it is possible to extract large-scale asymptotic behavior of the correlation 
functions after an appropriate renormalization procedure which is needed to remove 
UV divergences. The functional formulation is advantageous also because the Green 
functions of the Fourier-decomposed stochastic MHD can be calculated by means of 
Feynman diagrammatic technique. 

Action (5) is given in a form convenient for a realization of the field theoretic 
perturbation analysis with the standard Feynman diagrammatic technique. Free 
(bare) propagators A can be found from the quadratic part of the action (5) written 
in the form — (l/2)$/C<l> and by using the definition K.A = 1, where 1 denotes the 
diagonal matrix whose diagonal elements are the transverse projectors (our fields are 
solenoidal). One obtains 

\ 

(8) 





AW 







A vv 

3 s 






A bb 

3 s 




A v ' v 
j s 













A b ' b 

3 s 






A bb ' 

3 s 





J 

with the elements (wave-number-frequency representation) 



A-(k )W ) = A J V(-k,-o;) = 
J J —ioj + v k 

Af; ( k )W )=A^ ( -k,. p - (k) 



— iu> + uq uq k 2 ' 

Attn \ 2 3 , 2 fffalO fc ~ 2 + 9b20 p 

The model has three triple (interaction) vertices 

-v'(v-V)v = v'jVjkiVkVi , (10) 

-t7'(6-V)6 =v' j V jkl b k b l , (11) 

6'[(6-V)t» - (v-V)b] = b^Ujuhvi , (12) 

where the tensor structure of the vertices in wave-number-frequency representation 
are 

Vjki = i(SjkPi + Sjipk) , Ujki = i(5jiPk - SjkPi), (13) 
where momentum p is flowing into the vertex via the auxiliary fields v' , and b' . 
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3. Renormalization 

3.1. Divergences of the model 

It can be shown [17] that for any fixed space dimension d > 2, the superficial 
UV divergences can exit only in the following one-particle irreducible (1PI) Green 
functions: T vv , T bb , and T v bb . They lead to local counterterms of the form oc i/V 2 t>, 
cx 6'V 2 b, and cx v'(b-\7)b which are already present in the action (5), therefore, the 
model is multiplicatively renormalizablc (the analytic terms in k 2 proportional to g v2 o , 
and g V 2o in (3), and (4) are not needed in this case, and the model can be formulated 
without them). 

The situation is more complicated in the two dimensional case, where additional 
UV divergences appear. They are related to the 1PI Green functions T v v , and T b b . 
In this situation the formulation of the model without local (analytic in k 2 ) terms 
cannot give, in general, multiplicatively renormalizable model because the nonlocal 
terms of the action is not renormalized since the divergences produced by the loop 
integrals of the diagrams are always local in space and time (see, e.g., [13]). Thus, the 
simplest way how to restore the renormalizability of the model (or how to include the 
corresponding local counterterms cx v'W 2 v' , and cx b'W 2 b' in the renormalization) is to 
add corresponding local terms to the force correlation functions. It is shown explicitly 
in (3), and (4). In language of classical hydrodynamics the forcing contribution cx k 2 
corresponds to the appearance of large eddies convected by small and active ones and 
it is represented by the local term of v'V 2 v'. In its analogy the term b'V 2 6' is added 
to the magnetic forcing. 

Thus, in two dimensions, the model (5) is renormalizable by the standard power- 
counting rules, and for limits e — > 0, 8 — > possesses the ultraviolet (UV) divergences 
which are present in five aforementioned 1PI Green functions. It means that model is 
regularized using a combination of analytic and dimensional regularization with the 
parameters e, and 8 = (d—2)/2. As a result, the UV divergences appear as poles in e, 
(5, and their following combinations: 2e + 8, and (a + l)e + S. The UV divergences may 
be removed by adding needed counterterms to the basic action Sb which is obtained 
from the unrcnormalized one (5) by the substitution of the renormalized parameters 
for the bare ones: 3^10 -» fi 2e g v i, 9via -> \iT 2b g v2 , gbw —> fJ? ae gbi, S&20 -» H~ 2S gb2, 
vo — *• v-, mo — > u, where \i is a scale-setting parameter having the same canonical 
dimension as the wave number. 

In what follows, we shall work with, in our case the most convenient, minimal 
subtraction (MS) scheme, i.e., we are interesting only in the singular (pole) parts of 
divergent 1PI Green functions which are included in the renormalization constants. 
They give rise to the counterterms added to the basic action to make the Green 
functions of the renormalized model UV finite. In our model, the counterterms have 
the form 



where the renormalization constants Z^, i = 1, 2, 4, 5 renormalizing the unrenormalizcd 
parameters cq — {g v io,gv20,gbiOj9b20i l/ a, u o}, an d the renormalization constant Z3 
renormalizing the fields b, and b'. They are chosen to cancel the UV divergences 



count 




+ \(Z 4 - l)uv 3 g v2 pT 2& v>X7 2 v> + \ (Z 5 - 1) u 2 v 3 g b2 ^ 2S b'\7 2 b> 
+ (1-Z 3 )«'(6-V)6] , 



(14) 
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appearing in the Grren functions constructed using the basic action. The remaining 
fields v' , and v are not renormalized due to the Galilean invariance of the model (5). 

Renormalized Green functions are expressed in terms of the renormalized para- 
meters 

9vi — 9vw f- 2e ZfZ 2 , g V 2 — 9v2o f 26 Z\ZiZ± 1 , 

9bi = 9bio ^ 2at ZxZ^Z^ 1 , g b2 — 91,20 V? 8 Z^Z^Z^Z^ 1 , (15) 

v = vq Z^[ l , u = uqZ^Zi 

appearing in the renormalized action S R connected with the action (5) by the 
relation of multiplicative renormalization: S R {e} = S{e a } . The renormalized action 
S R , which depends on the renormalized parameters e(/i), yields renormalized Green 
functions without UV divergences. The RG is mainly concerned with the prediction 
of the asymptotic behavior of correlation functions expressed in terms of anomalous 
dimensions 7^ by the use of (3 functions, both defined via differential relations 



1i =i" 



d In Zj 



dg 



These definitions with expressions (15) yield the 7-functions 



with g = {g v i,g V 2,9bi,9b2,u} . (16) 



-f gvl = -271 


- 72 , 


7s61 = 


-71 - 2 72 + 73) , 


1gv2 = -271 


- 72 + 74 , 


7 3 62 = 


-71 - 2 72 + 73 + 75 


lv = 71 , 


1 

76 = ^73 , 


7u = 


-71 + 72) 



(17) 



and then /3-functions 

0gvi = 9vi (~2e + 271 + 72) , Pgbi = 9bi (-2 ae + 71 + 272 - 73) , 

Pgv2 = 9v2 (2(5 + 27i + 72 - 74) , Pgb2 = 9b2 (2<5 + 71 + 272 - 73 - 75) (18) 
0u = u(7i - 72)- 

3.2. RG equations 

Correlation functions of the fields are expressed in terms of scaling functions of the 
variable s = (k/fj,), s e (0,1). Then the asymptotic behaviour and the universality 
of MHD statistics stem from the existence of a stable IR fixed point. The continuous 
RG transformation is an operation linking a sequence of invariant parameters g~(s) 
determined by the Gell-Mann-Low equations 

TT-^ = (V( s )) with the abbrevation g = {g vl ,g v2 ,9bi,9b2,u} , (19) 
a Ins 

where the scaling variable s parameterizes RG flow with the initial conditions g\ s =i = g 
(the critical behaviour corresponds to IR limit s — ► 0). The expression of the (3(g(s)) 
function is known in the framework of the 5, e expansion (see (24) and also (18)). The 
fixed point g*(s — > 0) satisfies a system of equations (3 g (g*) = 0, while a IR stable 
fixed point, weakly dependent on initial conditions, is defined by positive definiteness 
of the real part of the matrix ft = (d(3 g /dg)\ g * (the matrix of the first derivatives 
taken at the fixed point). In other words, a fixed point is stable if all the trajectories 
g(s) in its vicinity approach the value of the fixed point. 

The initial conditions <7| s ->i = g of the equations (19), dictated by a micromodcl, 
are insufficient since our aim is the large-scale limit of statistical theory, where 
g* = g\ s ^o- As was mentioned already, the RG fixed point is defined by the equation 

/?(<?*) = 0- (20) 
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For g(s) close to g* we obtain a system of linearized equations 

(j*^-n) (g-g*) = 0, (21) 

where I is (5 x 5) unit matrix. Solutions of this system behave like ~g = g* + 0(s^ j ) 
if s — > 0. The exponents £j are the elements of the diagonalized matrix f2 dm9 = 
(£i> C27 £3, £4>£s) an( i can be obtained as roots of the characteristic polynomial 
Det(f2 — £7). The positive dcfineteness of represented by the conditions Rej(£) > 
0, j = 1, 2, ...5 is the test of the IR asymptotical stability of discussed theory. 



3.3. One-loop order calculation 

In the standard MS scheme [35] the renormalization constants have the general form 

Zi = 1 - FiP S '\ (22) 

where the terms P S ' e are given by the linear combinations of the poles and the 
amplitudes are some functions of g v \, g V 2, <7bi, 562; and u, but are independent of 5 
and e. The amplitudes Fj = F^F^ are a product of two multipliers F^\F^ 2 \ One 

of them, say, F^ is a multiplier originating from the divergent part of the Feynman 

(2) 

diagrams, and the second one, F} is connected only with the tensor nature of the 
diagrams (see discussion in [29] for details). 

It can be explained by the following simple example [29] (the example is taken 
from a problem with anisotropy, i.e., where another arbitrary unit vector n exists but 
the conclusions are the same). Consider a UV-divergent integral 



I ( K .11 = n.UjKiKm / <Z q + ^2)1+25 ^ ~i 



(summations over repeated indices are implied) where m is an infrared mass. It can 
be simplified in the following way: 

r°° 2 q zs 
7(k,n) = niUjkikmSijim dg 2 ( g 2 + m 2)i+2<5 > 



where 



Sijlm = ^ ^_ ^ (Sijdlm + SilSjm + 5 im 5jl — -{SijSi rn + SuSj m + Si m Sjl)), 



f 

Jo 



dq 2 i 2S r(5 + i)r(«y) 



/„ * 2( g 2 + m 2 ) 1 + 2 < 5 2m 25 r(2(5+ 1)' 

and Set — 2n d / 2 /T(d/2) is the surface of unit the d-dimensional sphere. The purely 
UV divergent part manifests itself as the pole in 25 = d — 2; therefore, we find 

UV div. part of I = —AF[ 2) k 2 + F, (2) (nk) 2 ), 
26 

where F[ 2) = F 2 (2) /2 = (1 ~ d)S d /M(d + 2) (F x (1) = F 2 (1) = 1). It has to be mentioned 
that in spite of the above simple example in our calculation we shall introduce the 
needed IR regularization by restriction on interval of integrations. 
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(2) (2) 

In the standard MS scheme one puts d = 2 in i 7 ^ , F 2 , therefore the d- 
dependence of these multipliers is ignored. As was discussed in [29], for the theories 
with vector fields and, consequently, with tensor diagrams, where the sign of values of 
fixed points and/or their stability depend on the dimension d, the procedure, which 

(2) (2) 

eliminates the dependence of multipliers of the type F± , F 2 y on d, is not completely 
correct because one is not able to control the stability of the fixed point when d = 3. 
Therefore, in [29] it was proposed to slightly modify the MS scheme in such a way 
to keep the d-dependence of F in rcnormalization constants Zi. Then the subsequent 
calculations of the RG functions (/?— functions and anomalous dimensions 7j) allow 
one to arrive at the results which are in qualitative agreement with the results obtained 
in the framework of the simple analytical regularization scheme, i.e., one is able to 
obtain the fixed point which is not stable for d = 2, but whose stability is restored for 
a borderline dimension 2 < d c < 3. In what follows, it will be shown that it is really 
our case, thus we shall apply this modified MS scheme in our calculations. 

Now we can return and continue with RG analysis. Using the RG routine the 
anomalous dimensions Jj(g v i, g V 2, 5&i, 9b2) can be extracted from one-loop diagrams. 
Thus, the extraction of the UV-divergent parts from one- loop diagrams gives Z- 
constants in the form 

Sd T , (9v2 9vi\ , , {9b2 gbi\ 



Z 1 = l + 



(2tt) c 



Zo = 1 



z 3 = i + 



Sd 



(2TT) d (u+ 1) 

Sd . f 9vi 



Ml 



9vi 
2e 



25 



9bi \ 
2ae) 



(27T) 

Sd 



9v2 

28 



9bi 9b2\ 
2ae 25 J 



(23) 



A 4 



z 5 = i + 



(2n) d g v2 \2S + 4e 
Sd A 2 



2ug vl g v2 
2e 

g v i gbi 



U 9 2 v2 

26 



9bi 



25 



+ 9vi gb2 



(2n) d (u+ 1)56 2 \25 + 2e(l + a) 2e 
and, in consequence, the lowest order 7-functions are 



f 4ae 

9v2 9bl 

2ae 



2gb\gb2 
2ae 

_ 9v2 gb2 

25 



9b2 

25 



7i 



Sd (uX 5 9v + A 6 5fc) 
73 = S d A 7 {-g v + gb) , 

~ A. a.. 

75 = s d 



72 = s d 



(Ai g v + A 3 g b ) 



1 



A 4 



74 = S d —{ugl+gl). 
gv2 



(24) 



9b 



(1 + u) g b2 



where Sd = 
9vi + 9v2, 9b 



Sd/(2n) d , Sd denote d-dimensional sphere Sd = 2n d / 2 /T(d/2), g v = 
= 961 + 362 , and A-coefficients depend only on dimension d: 



Ai 



A fi = 



d-1 
2d ' 
d-1 

4(d + 2) 



A 2 



A fi = 



d-2 

2d ' 
d 2 +d 



A, 



d-3 
2d 



A 4 



Ad{d + 2) ' 



4d(d + 2) 



A 7 = 



d(d + 2) 



(25) 



Substituting (24) into /3-functions (18) one can obtains /3-functions in the one-loop 
order approximation. Note that in two dimensions the 7-functions are 

I i2\ 1 (g v - gb) m I 



(2) 
7i 



32tt 



(ug v +g b ) , 



(2) 
7 2 



8tt (m+1) 



^j^(-9v+g b ), 
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7 f> = -LM±J?!), 7 ( 2 ) =0 (26) 

OZ7T 3^2 

and, in correspondence with [27] Z5 = 1, which is a specific property of the two- 
dimensional MHD turbulence because there are no UV divergences in the 1PI Green's 
function T b b in the one-loop approximation. Here we emphasize that in general case 
of d dimensions 75 7^ and Z5 7^ 1 . 



4. Fixed points 

Case of passive vector admixture 

Here we briefly consider the case when the magnetic field can be treated as a passive 
vector field in the developed HD turbulence. Notation the "passive" magnetic field 
means that the Lorentz force acting on conductive fluid can be neglected at large 
spatial scales, thus, the Lorentzian term (b • V)b in the Navier-Stokes equation can 
be omitted. Just then the vertex function T v bb is finite and the term containing 
Z3 in S coun t does not exists. Therefore, the magnetic field is not renormalized and 
73 = 0. Furthermore, some diagrams of T v " , T v v ,T b b containing the vertex T v bb can 
be omitted and Z-constants as well as 7-functions are reduced. Resulting 7-functions 
take the form 

9v 



71 = S d u A 5 g v , 72 = S d Ai 



u+ r 



74 = SdA 4 — g v , 75 = S d7 — — T . (27) 

9v2 562 

Substituting of 7-functions (27) and 73 = into /3-equations (18) one obtains a system 
of fourth nonlinear equations (3 gv \ — /3 gv2 — (3 g bi = P g b2 = for 3, and one equation 
p u = for u. The last one gives u* = 0, or, nonzero universal inverse Prandtl number, 




(28) 

In the first case of u* = one obtains only two fixed points (with zeroth 3^, <7^ 2 ) : 

1. g*i — 0,3*2 — — 25/XiSd which is non-physical (negative), and, 

2. g*! = 2e/AiSd, <7* 2 = which in unstable. 

Let u is given by (28). Then apart from the Gaussian fixed point g*^ = g* 2 = 
9bi = 9b2 = 0' with no fluctuation effect on the large-scale asymptotics, there are 
following fixed points with g£ 2 = 0: 

_ 2(rf-2KV + l) ~-i _ 
U ' 9v2 ~ 2d 2 -3d+2 bd ' 5bl ~ U ' 
4ed(u* + l) ~-i 

3(^_^ Arf ' 9v2 - U, fffti - 0, 

4e(3d 3 + d 2 (4e - 9) - 6d(e - 1) + 4e)(u* + 1) 
9(d+2e-2)(d- l) 2 
. _ 86 2 (d 2 -2)( M * + l) ~-i , 

9(d+2e-2)(d-l) 2 ^ ' 561 _0 - 

Next three fixed points are the same as the last (l*)-(3*) with different 3j* 2 : 
(la*) g* b2 = (d 2 -2)/d(d-2): 



(1*) 


9vi 


(2*) 


# 

9vi 


(3*) 


* 
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(2a*) ee (3a*) g* h2 = 3(d - l)(d + 2e - 2/2(d - 2)e. 

The points (2a*) and (3a*) have the same because g* t of the point (2*) is equal 
to the sum (g* x + g* 2 ) of the point (3*). Note that g* 2 has discontinuity at d = 2. 

The "thermal" point (1*) is generated by short-range correlations of the random 
force [27] and has negative g* 2 . The second fixed point (2*) is unstable. The physical 
meaning has the third "kinetic" point (3*) whose parameters {gi,g2,u} dependence 
on the dimension d is shown in figure 1. for physical value of e = 2. 




Figure 1. Dependence of the parameters {g v i, 9i>2,«} 

on the dimension d for e = 2 at the kinetic fixed point (29). 



Setting e = 2 and u* from (28) one obtains 

* _ (27r) d 8( M * + l)(3d 3 -d 2 -6ci + 8) , _ (27r) d 32(u* + l)(d 2 - 2) 
9vl S d 9(d-l) 2 (d + 2) ' 9v2 S d 9(d-l) 2 (d + 2) ' [ ' 

In this case the sum of 5* x +3*2 = .9* = (27r) d 8d(u* + l)/3(d— 1)5^. Detailed numerical 
calculations have shown that the region of stability of this point is limited by the value 
of parameter a < 1 and this limiting value does not depend on the dimension d. 
This stable region is denoted as region A in figure 2. 



4-2. Case of active vector admixture 

In the full self-consistent system, the RG equations yield besides the known fixed point 
in the kinetic regime also the nontrivial magnetic fixed point. If the both are stable in 
the same region of parameters then the choice between two possible critical regimes 
will depend on initial conditions for RG equations, i.e. critical behavior of the system 
is non universal. 



4-2.1. Kinetic fixed point. The nontrivial stable kinetic fixed point of RG equations 
has been found to be the same as in the previous case of passive magnetic field 
admixture because the (3- functions f3 gv i,f3 gv 2 are the same for zero gti,gb2- Only 
difference was found in the stability region in dependence on parameter a: the stable 
region is enlarged by new region B unlike the case of passive magnetic field admixture, 
see figure 2. The critical dimension d c continuously decreases from 3 to 2 in dependence 
on value of parameter a from the interval (1.1595, 1.427). It confirms the results of 
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Figure 2. Stability regions of kinetic point and the critical dimension d c 
dependence on parameter a. The region A spreads down to a = 0. 

[27] that the stability of kinetic scaling regime is strongly affected by the behavior of 
magnetic fluctuations. 

Figure 3 shows values of the charges g v i,g V 2 which continuously depend on value 
of nonzero e < 2, for two special case of d equal to 2 and 3. The right axle corresponds 
to the physical value of e = 2. While the both charges remain nonzero (positive) for 
d = 2, in three dimensions one of them, g V 2 rapidly decreases for e — > 0. The stable as 
well as unstable regions depends on parameter a and the critical value of a remains 
the same for e = 2 following from figure 2, or greater for e < 2 (the critical a increases 
for e -» 0). 




Figure 3. Dependence of the parameters {g v \ , g v 2, u} on value of e 
for 2- and 3-dimensions at the kinetic fixed point in the general case. 
Dashed line shows the critical value of a at the stability region limit. 



4-2.2. Magnetic fixed point. We have shown in (26) that in two dimensions the 
function 75 vanishes and then both functions (3 g bi and I3 g b2 contain the same linear 
combination of 7 functions. Thus, at least one of the magnetic charges (ffbi, ^2) must 
be zero in fixed point. But in the other dimensions this restriction does not take place. 

Here we restrict ourselves only by finding nontrivial magnetic fixed point. In 
[27] there was mentioned that it is characterized by zero g* x and u* . Therefore, the 
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set of five equations of zero /3-functions (18) is reduced to three equations. Applying 
g v \ = u = in (18), (24) and (25) one obtains the set 

aig v2 + a 2 g 2 v2 + a 3 g v2 gb - a 4 g b = , 
-Ao + a 5 g v2 + a e g b = , 

aifita! + a 5 g v2 gb2 + a 6 g b2 g b - a 7 g v2 g b = , (30) 

where 

2ae 2(d-2) (d - 1) 

_ (d 2 - 5) (d 2 - 2) _ (d 2 + d - 1) 

fl3 ~d(d + 2)' a4 ~4d(d + 2)' a5_ d(d + 2) ' 

(5d 2 -3d-32) (d-2) 
fl6= 4d(d + 2) ' ^HsT" (31) 

Positive coefficients 01,07 vanishes at d = 2, 03 and 06 are positive for d > 2.236 
and d > 2.848, respectively. The set (30) can be analytically solved with respect to 



2.8 
2.4 



2.0 




y/ 1 

(0.723, 3) / 




unstable 


(0 




stable 


on-physi 


(0 - 146 ' 2 " 362) (1.2a/ 
unstable 


c 



0.4 



0.8 a 



1.2 



Figure 4. The stability region of the magnetic fixed point 
in the plane of {d, a} for the physical value of e = 2. 

g V 2,gbi,gb2- Because all g% must be positive, the system (30) with g v \ 
the only solution, 



u = gives 



g V 2 = 



where 



9b2 

9b = 
D = 



h ~ a&9b 
a 5 

gb{a 6 gb - A ) 



561 = 9b - 



o 5 a 7 



9b{af>9b - A ) " 0507 
05(01 + 2a & g b - A ) 



05(01 + 2o 6 56 - A ) 

-oio 5 a 6 + o 3 a 5 Ao - 2a 2 a e A + a 5 \/D 

2(a4a| + 03a 5 a 6 - o 2 a^) 
a\a\ + 4010405^0 + 2oia 3 o 6 ^4o + a 3^1o + 4o 2 04^4g . 



(32) 



(33) 



Note that the parameters a and e appears in the solution only as the product at 
in Ao. The physical value is restricted by the inequality oe < 2. Denominator in 
the expression (33) for g b is zero at do = 2.2628 (and it is positive for d > do), 
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therefore, at d we can expect discontinuity and/or divergence. Numerical analysis 
of the expressions (32) shows that all gi have a discontinuity at do, and, a physical 
solution cannot exist for any a, e if d < d . The stability region of the magnetic fixed 
point and the corresponding critical dimension d c was determined numerically and it 
is shown in figure 4. 




Figure 5. Dependence of the parameters {g v 2, gbi, 9m} on the dimension d 
for e = 2 at the magnetic fixed point (32) for a = 1, 0.5, 0. All have 
discontinuity at d = 2.2628 (chained vertical line). 

Figure 5 demonstrates the charges g V 2,gbi,9b2 dependence on dimension d for several 
values of parameter a. First, we have found that g V 2,gb2 tend to infinity at limit 
value do- For increasing dimension d from 2 up to d the charge gb2 increase from a 
small positive value up to infinity at do and, therefore, g v2 decrease here from a small 
negative value to minus infinity at do (because g V 2 oc —a e gb and both and gb are 
negative in this dimensions). The charge gbi rapidly decrease to zero at d — 2.352 
for decreasing d and continue to minus infinity at do- These limiting value are in 
correspondence with numerical calculation of the stability region - the system losses 
stability for the critical dimension d c lower than approximately 2.36 for arbitrary 
parameter a. 
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5. Discussion and conclusions 

In this paper we revised the calculations of stability ranges of developed 
magnetohydrodynamic turbulence in the frame of double expansion scheme. The 
modified standard minimal subtraction scheme [29] has been used in the dimension 
region of d > 2 up to d = 3 in both cases of the magnetic field treated as a passive 
as well as active vector admixture. We confirm existence of the known " kinetic" fixed 
point (corresponding to Kolmogorov scaling regime) what is the same in the both 
considered cases and only difference is in the stability region: the critical dimension 
d c is achieved for a slightly higher value of a-parameter of a magnetic forcing in the 
case of active magnetic field. Limit value of the inverse Prandtl number at d = 3 
restores the value of u = 1.393 which is known from usual e-expansion, and it fluently 
rises to u = 1.562 at d = 2, (figure 1). 

It was believed earlier that in the double expansion being defined for the space 
dimension to be closed to two the results obtained in two dimensions can not be 
applicable to opposite dimension interval end closed to three. Here we have showed 
that the double expansion in exact d-dimcnsional formulation gives some critical 
dimension d c above which the scaling regime is governed by the competition of the 
stable kinetic and magnetic fixed points which exists in three dimensions. 

A new nontrivial results of the present paper is connected with derivation of 
the exact analytical expression for the nontrivial "magnetic" stable fixed point with 
u = g v \ =0 but nonzero g V 2,9bi and gm as well as specification of the borderline 
dimension d c . A physical region of the RG fixed point lies below the ae = 2 line, 
see in figure 4. This point completely losses stability below the critical value of 
dimension d c — 2.36 (independently on the a-parameter) and also below the value 
of a c — 0.146 (independently on the dimension). Thus we confirm, in particular, that 
thermal fluctuations of the magnetic scaling regime may occur, and, in comparison 
with earlier results our value of the borderline dimension (d c — 2.36) is significantly 
lower than in the e expansion [10] (d c = 2.85) and rather lower than in the 'modified' 
double expansion introduced in [27] (d c = 2.46) but it is rather higher then value 
(d c = 2.2) calculated in the frame of the McComb's renormalization [21]. 

Note that the stability of any regime determines the concrete Alfven ratio ta 
(ratio of kinetic and magnetic energy density in MHD turbulence, see [22, 23], for 
example). Once the stationary scaling regime becomes and stands, the Alfven ratio 
is fixed (i.e., it means that the fixed point is reached in the field RG terminology). 
Thus the injected energy necessary to steady the stationary scaling regime must have 
specific value, or, in another words, all "coupling constants" gi are fixed in scaling 
regime with values which arc dependent on dimension d. In like manner the inverse 
Prandtl number u = njv {r\ is magnetic resistivity) is thus fixed. Vcrma [23] has 
obtained r\jv = 0.85/0.36 = 2.36 in 3-dimensions for large ta ~ 5000 (corresponding 
to region of the kinetic regime) and for zeroth normalized cross-helicity. For smaller 
ta this ratio decrease to 0.69 for ta = 1, and, the both rj and v vary approximately 
as d~ x l 2 [22]. We have mentioned above that in our double expansion calculation in 
the kinetic point we have fixed the ratio u = rj/v with its (/-dependence showed in 
figure 1. The magnetic fixed point is characterized by decreasing value of u to zero 
what is in correspondence with results of [23]: his calculation gives for decreasing ta 
(magnetic regime) in 3-dimensions also decreasing value of rj/v as one can expect in 
the magnetic fixed point. 
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